function [L y ok]=putere_directa(A,y,max,tol)
%	Functia care calculeaza valoarea proprie cea mai mare a matricii A, precum si
%vectorul propriu corespunzator acestei valori.
%
%	Date de intrare:
%		-A -> matricea a carei valoare proprie se cere de determinat;
%		-y -> vectorul propriu initial;
%		-max -> numarul maxim de pasi;
%		-tol -> toleranta, adica || A*y - L*y || (L->valoarea proprie)
%
%	Date de iesire:
%		-L -> valoarea proprie determinata;
%		-y -> vectorul propriu determinat;
%		-ok -> o valoare care este 1, daca || A*y - L*y || <= tol si 0 altfel

n = 0;

L = 0;


while( n <= max && (norm(A*y - L * y) > tol) )
	n = n + 1;
	z = A*y;
	y = z./norm(z, 2);
	L = y'*A*y;
endwhile

ok = (norm(A*y - L * y) < tol);

endfunction
